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Noise induced jumping between meta-stable states in a potential depends on the structure of the 
OS ' noise. For an a-stable noise, jumping triggered by single extreme events contributes to the transition 

as : probability. This is also called Levy flights and might be of importance in triggering sudden changes 

Os . in geophysical flow and perhaps even climatic changes. The steady state statistics is also influenced 

by the noise structure leading to a non-Gibbs distribution for an a-stable noise. 
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I. INTRODUCTION 



Noise induced jumping between meta-stable states separated by potential barriers is common in physical systems. 
The time scale for the barrier penetration depends on the structure of the noise. Most often the noise is Gaussian. 
However, non-Gaussian noises distributed with power-function tails, Levy flights, are observed in many different 
physical systems |lj such as turbulent diffusion |^,3| and vortex dynamics [Q. Levy flights also seems to be a common 
feature in dynamical models || an d critical phenomena || . 
. . The Levy flights can result from a Langevin equation driven by a-stable noise and give rise to anomalous diffusion 
of a random walker with position r(t) such that (\r(t) ~ r(0)| 2 ) cx Dt 2 l a where D is a constant and < a < 2 Jl6| . 
JL| The case a = 2 corresponds to normal diffusion where D is the diffusion constant. The exponent a is related to the 
scaling of the tail of the probability distribution for the increments of the random walker, P(X > r) oc r~ Q . For a > 2 
!— 1 1 the second moment exists and by the central limit theorem the random walker reduces in the continuum limit to a 
Gaussian random walker unless the diffusion takes place on a fractal set like in a quenched random medium j7j. In 
this case the random walk can be sub-diffusive. Another example of a process which can be sub-diffusive is the Levy 
^v^j . walk where a random walker has a constant speed in between discrete stochastic time points (a renewal process) with 
a power-function tail distribution. Note that since the time process is discrete for a Levy walk it cannot result from 
f^) , a Langevin equation. 

t-H . Anomalous diffusion was first observed in hydrological time-series ||. Recently evidence for a-stable statistics in 
atmospheric circulation data has been reported [fl0|| . In a long paleoclimatic time-series an a-stable noise induced 
jumping in a double- well potential was found Q. In both cases a was found to be around 1.7. The latter describes 
a jumping, in glacial times, between two climatic states governed by the oceanic flow forced by random fluctuations 



from the atmosphere. Understanding the role of extreme events and the time-scales for these climatic shifts is the 



main motivation for this study. 



In this paper we will interchangeably use the physics jargon, (x), and the mathematics jargon, E[x], for the 

j^n *- expectation value for x. The latter will be used in the case of conditional expectations. We use the usual convention 

that probability distribution functions, P, are capitalized and probability density functions, p — dP/ 'dx, are in small 

• • ■ letters. 
> ■ 

jj ■ II. THE q-STABLE DISTRIBUTIONS 

d; 

For distributions with power-function tails, P(X > x) oc a;" 7 , only moments of order less than 7 exists = 
00 for/3 > 7). For < 7 < 2 a generalized version of the central limit theorem applies, namely that the average of n 
independent stochastic variables from the distribution P asymptotically will have an a-stable distribution as n — > 00 
with a — 7. The a-stable distributions are defined by their characteristic functions, (exp(ikX)) = exp(— a a \k\ a /a). 
The a-stable distributions are stable with respect to averaging, Y n = n~ x l a X)"=i ^ii meaning that Y n has the same 
distribution as Xi where the Xj's are i.i.d. (independent identically distributed) a-stable, thus the phrase 'a-stable'. 
As for the case of Gaussian noise, the dynamics of a noise driven system with power-function tail distributions for 
the noise increments, P{X > x) oc x~",0 < a < 2, will reduce to a system with an a-stable noise in the continuum 
limit, described by a Langevin equation B, 
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dX = f(X)dt + a{X)dL, 



(2.1) 



A random walker with a-stable noise increments will be super-diffusive due to the large jumps from the tails of the 
distribution surviving the averaging in the continuum limit. See Appendix A for a further short description. 



III. THE FOKKER-PLANCK EQUATION 



The probability density for X in (2.1) is determined from the Fokker-Planck equation (FPE), see Appendix B for 
a derivation, 



d t p(x) = -d x [f(x)p(x)] 
' lkx ^(k - k 1 )\k\ a p{k 1 )dkdki. (3.1) 



The second term on the right hand side is expressed in terms of the Fourier transformed probability density, p{k). 
This term reduces to the ordinary diffusion term, d x 2 [a 2 (x)p{x)]/2. when a = 2. In this case the solution for the 
stationary probability density function can be expressed explicitly in the well-known form, 

P (,)c<^ex p{ 2^^}, (3.2) 



For a < 2 the FPE (3.1) is non-local in spectral space. This is a reflection of the super-diffusivity of the process 
( [2.l|) . Besides the Gaussian case we can only solve the FPE explicitly for a = 1. This is the case of a system driven 
by Cauchy distributed noise having the probability density, q(x) — l/[7r(l + a; 2 )], se e A ppendix C for further details 



on the Cauchy distribution. We are using p(x) for the probability density for X in (2.1) and q(x) for the probability 
density of the noise. Then the stationary FPE becomes, 

i J f(ki — k)p{k\)dk\ = sgn(k)ap(k) (3.3) 

where the noise intensity a is taken to be constant. From taking the derivative with respect to k on both sides of 
(|3.3[) and performing a partial integration on the left-hand side it follows that 



i J f(h - k^^ik^dh = sgn(fc)ap< m )(/c) (3.4) 
for any m. The solution is p(k) = e~^' k ', where A is determined by 

i J /(fci - k)e- x{kl ~ k Uk x = if{i\) = sgn(fc)cr. (3.5) 

Thus the solution is determined by the analytic continuation of f(x) into the complex plane, provided it exists. Note 



that the solution also apply for k = where the r.h.s of (3.4) jumps, since from the def inition of the Fourier transform 



of the probability density we have p(0) = E[l] = 1. By com plex conjugation of (3J3) we get if{—i\*) = sgn(— k)a, 



so for k < the solution is given as —A*, where A solves (3.5) for k > 0. With A = (3 + iS the characteristic function 
is given as pik) = e -P\ k \ e - lSk ^ p or ^(fc) to be a characteristic function we must have ft > 0, and the stationary 
distribution is 



^ 1 0i 



fr[ n/3f + (x + 5i) 2 

where the sum is over the N zero points of the complex function if(i\) — sgn(fc)cr in the upper half-plane {(3 > 0). In 
this solution of th e stationary FPE there is an indeterminacy since any p{x) with ^ . pi = 1 is a probability density 
that satisfies (|3.l|) . 

The indeterminacy might be related to the problem of conservation of probability. If there is a finite probability 
for the random walker to escape to infinity, it must be reinserted into the system for a stationary probability density 
t o b e conserved. Then the indeterminacy in the reinsertion could result in the indeterminacy in the coefficients pi in 
(|3.6[). However, the indeterminacy can be lifted in the limit a — > 0. When the intensity of the noise becomes small the 
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Cauchy distribution approaches a (^-distribution (when acting on functions that are bounded by \x\@ for some (3 < a 
as x — > oo ). T hen we can approximate the system by a system with discrete states and the stationary Fokker-Planck 
equation ( [3.3] ) is approximated by N transition (Master) equations for the weights Pi,i = 1, ...N, 

Pi = ^jPjP(j -> 0) ( 3 - 7 ) 



where i,j represents the N minima defined in (3.6). The transition probabilities p(i — > j) are related to the transition 



waiting times which will be defined in the following. 



IV. THE POTENTIAL 



Before proceeding we will define the drift term as resulting from a potential. The governing equation then describes 
a massless, viscous particle moving in a potential, f(x) = —dU/dx. As an example for study we define the potential 

as 

17(a) = 4(a;/A) 4 + h{x/A) 3 - 8(a;/A) 2 - 3h(x/A). (4.1) 

U(x) is a double-well potential for —16/3 < h < 16/3. Ah is the level difference between two potential minima at 
x = —A = a and x = A = c. The local potential maximum between the two minima is at x = — 3/iA/16 = b, and the 
potential values are [U(a), [7(6), [7(c)] = [-4(1 - h/2), (3/i/16) 2 (8 - 3/i 2 /64), -4(1 + h/2)}. See figure 1. The results 
are readily generalized to other forms of the potential U(x). 
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FIG. 1. The potential ([O]). ill = 17(a) — U(c) is the potential difference between the two minima, w = b — a is the 'left 
half-width'. Units are arbitrary. 



V. WAITING TIME 



The waiting time for jumping between the two potential minima (from a to c) of U (x) defined above is exponentially 
distributed. With p ac (T > t) being the probability of staying in minimum a longer than t we have p ac (T > t) — 
exp(— t/T ac ) with a mean waiting time T ac . This follows from the Markov property of the Langevin equation in the 
discrete state limit, since we have 

P{t<r<t + At)/ At - (1 - X ac Atf At \ ac 

— > A ac exp(-A ac £), (5.1) 

as At — > 0, where A ac = l/T ac is the transition probability intensity. In the non-discrete case, a little more rigorous 
treatment is needed |T^ ]. However, the result holds, if the potential wells are substituted for the minima, and the 
waiting time is defined as the time between consecutive crossings of a and c. 
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A. Gaussian noise and Arrhenius formula 



In the case of Gaussian noise in (2.1) T ac can be calculated from the backward Fokker-Planck equation [521] , 



J-oo J a 



(5.2) 



and correspondingly for T ca . By using the saddle-point approximation on (5.2) we obtain the Arrhenius formula, 

T QC ocexp(2[C/(6)-f/(a)]/a 2 ). (5.3) 

For comparison with the case of a-stable noise figure 2 displays the standard result of a numerical simulation in 
the case of Gaussian noise. Figure 2 (a) shows the simulated process with the potential in figure 1. 2 (b) shows 
the simulated probability density function and the right-hand side of (|]^). Figure 2 (c) shows the time-scale for 
jumping as a function of the parameter h. The time-scale is calculated from the exponential distribution of times 
between consecutive crossings of the levels a and c. Figure 3 shows the number of crossings (from a to c and from 
c to a respectively) with a waiting-time larger than each waiting-time measured, normalized by the total number of 
crossings. These points are situated on straight lines in the semi-logarithmic plot where T ac and T ca are the slopes 
of the lines. Figure 2 (c) shows the time-scales for seven simulations with different h. The curves are the time-scales 
calculated from (|5.2|). 
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FIG. 2. A simulation of (2.1) with Gaussian noise and the potential shown in figure 1. (a) shows a realization and (b) the 



probability density function. The actual simulation is 1000 times longer than what is shown in (a). The smooth curve in (b) is 
the pdf calculated from ([j.2|). (c) shows the mean waiting ti mes, T ac (diamons) and T ca (triangles), for seven simulations with 
varying h. The curves are the waiting times calculated from (5.2). Units are arbitrary. 
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FIG. 3. The probability for waiting longer than t before jumping to the other well as a function of t obtained from the 
simulation. The slope of the upper curve gives T ca and the slope of the lower curve gives T ac - Units are arbitrary. 



B. a-stable noise 



In the case a < 2 the situation is radically different. The sample curves of the process are no longer continuous 
and the finite jumps or extreme events will contribute to the probability of jumping between the potential wells. The 
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probability, (X ac At + o(At)), for jumping from the left well, x < b to any y > b in a single jump in a time interva l A t 
is governed by the tail of the distribution, p(x) oc (x / a)^^ a+is> At/a. This is seen by observing that the process (2.1) 
can be obtained from the discrete process, X(t + At) = X(t) + f{X{t))At+[aAt 1 / a ]r ] {t) 1 for At -> 0, where rj(t) has 
an a-stable distribution with unit intensity. Thus we have 



X ac At re P(X(t + At) > b\X(t) < b)/P(X(t) < 6) oc 



p(u)du]p(x)dx 



b-a 



6-0 



p{u)du^ [{b-a)/<r]- a At 



(5.4) 



The inner integral is the probability of jumping from x < b to any y > b, and p(x) is the stationary probability 
density. The outer integral is dominated by the central part of the probability distribution. This result is exact in 
the At — > 0, a — > limit where p{x) — > 5{x — a). Thus we have 



T ac = c(a)[(b-a)/a] 



(5.5) 



where c(a) is some constant. So in this case we see that the waiting time scales with the 'left half-width' of the 
barrier, b — a = w, to the power a. The height of the barrier has no influence on the transition probability. The 
results are confirmed by numerical simulation. Figure 4 displays the numerical simulation using Cauchy noise, a = 1, 
and the same potential as in the case displayed in figure 2. Note the linear scale in figure 4 (c) showing the scaling of 
the time-scale with w. 
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FIG. 4. The same as figure 2 but with Cauchy noise, a = 1. Note the linear axis in (c). The curves are obtained from (5.5) 



VI. THE STATIONARY DISTRIBUTION 



Fo r the Cauchy noise driven system the indeterminacy in (|3.6[) can now be resolved by use of the Master equation 
(3.7). In the limit a — > the system can be approximated as a discrete two state system, with the two states 
corresponding to the two potential minima, at a and c. In this limit the system fulfill the stationary Master equation, 

=p a p(a —> c) -p c p(c — > a). (6.1) 

The transition probabilities are now p(a — > c) oc 1/T ac oc (b — a)/a and p(c — > a) oc (c — b)/a and we get, 

p a = l-Pc = (b-a)/(c-a). (6.2) 

Note that this is independent of exp(— 2[U(a) — C/(c)]/<t 2 ), which in the Gaussian case corresponds to the Gibbs 
distribution. Figure 5 shows the distribution p a , which is different from the Gibbs distribution, as a function of w . 
Figu re 4 (b) shows the probability density function from the simulation over-plotted the one calculated from ( |3.6| ) 
and (IO). 
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FIG. 5. The probability, p a = 1 — p c , for finding the particle in the left well as a function of w in the simulation with Cauchy 
noise. The curve is obtained using (S.2), The distribution is deviating strongly from a Gibbs distribution. 



VII. BARRIER PENETRATION 

When a is close to 2 we should expect the 'single jump penetration' of the barrier to become more and more unlikely 
and the continuous penetration dominating. The Levy decomposition theorem Q states that the cv-stable process 
can be decomposed in a Brownian process and a compound Poisson process. The 'continuous' barrier penetration can 
be estimated by considering the distribution to be truncated so that there are no jumps larger than the half-width 
of the barrier, w. The truncated probability for the noise, p t (x), is then defined by p t (x) oc p(x) for |x| < w and 
p t (x) = for \x\ > w. This part of the noise now has finite second order moment and we can estimate the variance as 
°eff oc J x 2 p(x)dx oc w 2 ~ a asymptotically for large w or small noise intensity a. The waiting time can be estimated 
as, 



T c ex a J exp(2[C/(6) - U{a)]/a 2 ^) 



(7.1) 



where c denotes 'continuous'. Note that this part of the process is not strictly continuous, since it contains jumps 
smaller than w. The time-scale for single jump penetration can be estimated from (5.5), 



Td c 
OC W 



(7.2) 



where d denotes 'discontinuous' and we have, 



^ « w 2 exp(-~c[U(b) - U(a)]w a - 2 ), 



(7.3) 



where c is a constant. So the relative importance of extremal jumping depends both on the heig ht and the width 
of the barrier. To illustrate the relative i mpo rtance of two jumping processes a simulation of (2.1) with an a-stable 
noise [ fl3[ with a = 1.7 and a potential (4.1) with h = 3, was performed. Figure 6 shows part of a realization of 
this process. Here it is seen that the jumping from the deep to the shallow well is governed by the discontinuous 
part, T d (c -> a) < T c {c -> a), while the jumping from the shallow to the deep well is dominated by the 'continuous' 
part, T d (a — > c) 3> T c (a — * c). Fo r pr oportioning the continuous and discontinuous processes in a given situation the 
prefactor and the constant, c, in (7.3) must be calculated or estimated. 
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FIG. 6. A realization of the process with a — 1.7. The potential used is shown in the insert. The jumping from the left 
(shallow) well to the right (deep) well is triggered by the - almost - normal diffusion. The jumping the other way is driven by 
the tail of the a-stable distribution, the extreme events. 

VIII. SUMMARY 

We have seen that the statistics of noise induced jumping between meta-stable states in a potential is different for 
a-stable noise from the usual Gaussian noise case. The stationary probability distribution deviates from the Gibbs 
distribution, and the waiting time for jumping depends in some cases more on the width than on the height of the 
barrier. This is the case where a single extreme event triggers the jumping. These observations might be of importance 
for understanding the triggering mechanisms of climatic changes, where the flow state of the ocean is trapped in a 
potential minimum, a stable climatic state. This flow is stochastically forced by the atmospheric flow. There are 
some evidence that this stochastic forcing is a-stable rather than Gaussian such that climatic shifts from one state to 
another could be triggered by single extreme events. This would perhaps explain why the climate models at present 
are not capable of reproducing the climatic changes observed in the geological records. The models are too coarse 
grained and contains to much diffusive smoothening to allow for extreme events. 
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APPENDIX A: THE ADDITION OF a-STABLE RANDOM VARIABLES. 

Textbooks on a-stable processes are now available jl4],[l3]], but for those readers not familiar with the a-stable 
distributions and processes a few notes are added in the following. 

When {Xi, i = 1, n} is a series of i.i.d. random variables, the distribution of the variable Y = c(n) Xj can 

be determined from the characteristic function, 

n 

{exp[ikY]) = (exp[ikc(n) Xj\) = 

(IT" eMikc{n)X 3 ]) = ( e *M»)*)« . (Al) 



If the distribution for Y is the same as for X the equation (Al) for the characteristic function, f(k) — (exp[ifcY]) is, 



f(k) = f(c(n)k) n (A2) 
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with the solution 



f(k)=exp[-a a \k\ a /a], (A3) 
c{n)=n- 1 / a . (A4) 

The c onst ant, a a /a is chosen so that it coinsides with the usual notation in the gaussian case, a = 2 . Only for a > 



does flAq) represent a characteristic function. It can be shown that the characteristic function (A3) corresponds to 
distributions with power-function tails, P(X > x) ~ x~ a For a > 2 the second moment of the distribution 

exists and sums of i.i.d. variables converges by the central limit theorem to the gaussian distribution, a = 2. For 
< a < 2 the distributions has a domain of attraction in the sense that sums of i.i.d. random variables with tail 
distributions, P(X > x) ~ x^ 1 ', under rather general conditions converges to an a-stable distribution with = 7. 
This is the generalization of the central limit theorem for a-stable distributions. The proof of this is similar to the 
proof of the central limit theorem for the normal distribution. It basically substitudes a limit, f(c(n)k) n — ► f(k) for 



(|A2|). The proof can be found in Fellers book, pp 574 - 581 |Lj 



Now we can intuitively understand the noise term, dL a , in the Langevin equation (2.1) as the continuum limit of 
addition of small increments, 

^ m 

AL a (At) = X(jAt/m) (A5) 

where X(t) is a random process with power-function tails, P(X(t) > x) ~ x~ a , and unit intensity. In the limit, 



m — > 00, AL a will be an a-stable noise. It follows from (A4) that dL a — dt 1 / 01 , which in the gaussian case is the 
well-known relation, dB 2 — dt. 

For a < 2 the a-stable variables have infinite variance. This concept can be difficult to comprehend when considering 
measurements from a given physical system. In the case a sample is taken, say of n measurements of the variable 
X, where X has an a-stable distribution with stability index a, then of course any of the measurements, x\, ■■■,x n 
of X is finite so that the sample variance, (x\ + ... + xfy/n, is some finite number. The variable Y = X 2 will 
have a tail distribution given by P(Y > x 2 ) — P(X > x) ~ x~ a = y~ a / 2 , so that, asymptotically for large n, 
Z n = n~ 2 / a (Yi + ... + Y n ) will have an a-stable distribution with stability index a/2. Imagine now that we estimate the 
(infinite) variance of variable X by taking samples of length n, estimating the variance as (Xf + ...+X 2 )/n = n 2 / a ~ 1 Z n . 
Then the estimate itself will be an a-stable process with stability index a/2 and intensity n 2 ^^ 1 . This estimate will 
be fluctuating with an intensity growing with n for a < 2. 



APPENDIX B: THE FOKKER-PLANCK EQUATION 

In the following the Fokker-Planck equation (3.1) corresponding to the Langevin equation (2.1) will be derived. 
The Fokker-Planck equation will be derived in spectral form using that the a-stable processes are defined by their 
characteristics functions. Following the lines of Stratonovich ]l7| we define the functional 

/= f R(y)d t p(x \y,t)dy = lim o I At (Bl) 

where 

lAt = At J R (yM x o\v,t + A t) -p( Xo \y,t)]dy (B2) 

R{y) is an arbitrary generator function, and p(xa\xi,t) is the conditional probability density at x\ corresponding 
to passing from xo to x\ during time t. Assuming stationarity we suppress the first temporal index, p(xo\xi,t) = 
p(xo, 0\xi, t) = p(xq, t\xi, T+t). For simplicity of writing we make the convention that J is to be read as (1/V2tt) ■ 
With p(xo\x,t) being a probability density in x we trivially have 



p(xo\x,t)dx = 1 (B3) 



and the Chapman-Kolmogorov equation 
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p(x \xi,t)= / p(x \x,T)p(x\x!,t- r)da 



(B4) 



For the functional we then get 



At 



I At = 

R (y)i I p(x\y,At)p(xo\x,t)dx-p(x \y,t)]dy 



At 



p(x \x,t){ / p(x\y,At)[R(y) ~ R(x)]dy}dx 



We now define the Fourier transforms 

R( x ) = f R(k)e lkx dk, R(k) = I R{x)e- lkx dx, 



(B5) 



(B6) 



p(x \x,t) = J p{x \k,t)e ikx dk, 
p(x \k,t) = p(x \x,t)e lkx dx 



similarly for f(x) in (2.1), and for a a {x) to be introduced below. However, for the probability density p(xo\x,t) we 
define 

(B7) 
(B8) 

consistent with the standard definition of characteristic function except for the factor v / 2~7r. With these definitions it 
is easy to derive the formula 

' f(x)p(x \x, ty kx dx = J f(h - kffixolk!, t)dh (B9) 
from which it directly follows that 

f(x)p(x \x, t) = J [J f(h - k)p(x \ki, ^dk^e-^dk 
Using the spectral representation for the generator function we get 

-ilAt - / P(x \x,t){ [ p(x\y,At)±- ( R{k)(e ik * - e lkx )dk}dy 



(BIO) 



^2vr 



At 



e ikx R(k)p(xo\x,t)^-E[e ik ^ x ^ +At ^ - l\X(t) = x]dkdx 



(Bll) 



The conditional expectation is evaluated using the Langevin equation (2.1) and the characteristic function of the of 
the alpha-stable Levy noise increment dL a , 



We get 



E{exp[ikdL a ]} = exp[-cr Q Ai|fc| Q /a] 



L£[ e *M*(*+ A *)-*] - l\X(t) = x] 



(B12) 



At 



_}_ E t e ik[f(x)At+o(At)+dL a ] _ ji 

At 1 1 

ik[f(x)At+o(At)]-a a At\k\ a /a _ -, 



At 



^ikf(x)-a a \k\ a /a 
as At — > 0. Substitution of ( B13| ) in (Bll) and combining with (B2) then gives 



(B13) 
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R{k)e ikx d t p(x \x,t)dxdk 
R(k)e lkx [ikf(x) - a a {x)\k\ a /a}p(x \x,t)dxdk 



(B14) 



Here we have permitted the scaling factor a a / a, corresponding to th e va riance a 2 dt of the noise increment in the case 
a = 2, to depend on the variable x. By eliminating the x by use of (B9) we get, suppressing xq, 



R{k)d t p{k,t)dk 



R(k)[ikf(h -k)- a a {k 1 - k)\k\ a / a)p{k u t)]dkdk x 



(B15) 



and finally since R(k) is arbitrary we get the spectral Fokker-Planck equation for the integrand, suppressing the t 
index, 



d t p{k) = / (*Jfe/(fci - k) - a a (h - fc)|/c| Q /a)p(fci)dfci 



Multiplying by e and using (BIO) gives the Fokker-Planck equation in the usual form 



d t p(x) = -d x [f(x)p(x)] - - 

a 



—ika 



a a (k- fci)|fc| Q p(fci)dfcdfci 



(B16) 



(B17) 



For the stationary Fokker-Planck equation the l.h.s of (B17) vanishes and the partial derivatives become total deriva- 
tives. The last term on the r.h.s is a generalized diffusion which formally can be written 



1 d a 

a dx a 



K(x)p(x)} 



a 



e~ ikx a a (k - kx^k^pik^dkdkx 



(B18) 



In the case a — 2 this is the usual diffusion term corresponding to Gaussian white noise excitation of intensity a 2 (x). 
For a < 2 the diffusion is non-local. The physical meaning of this term is that for a-stable processes there will, due 
to the fat tails of the distributions, be finite size jumps in the process. 



APPENDIX C: THE CAUCHY DISTRIBUTION 



The probability density can only be expressed explicitely for a = 1, and a = 1/2. For a = 1, Cauchy noise, the 
characteristic function is c{k) = exp(— <r|fc|) and its Fourier transform is, 

= — n n i \2i ( C1 ) 



For this distribution even the mean does not exist. Note that even though the density distribution (|Clj) is symmetric, 
p{— x) = p(x), this does not imply that (x) = 0. For a data sampling this manifests itself in the fact that the average 
of n data points, Z n = (X\ + ... + X n )/n is Cauchy distributed with the same intensity as Xj, so that there is no 
convergence for the series Z n , n= 1,2,..., it fluctuates exactly as the data Xi itself. 

A classical example of this characteristic of the Cauchy distribution is seen by considering the distribution of light 
on a line, L, from a point source, see figure 7. Since the light is uniformly distributed over the angles, p(9) — 1/tt, e 
[0, 7r], the distribution on the line will be, p(x) = p(0){d0 / dx) — l/[<$7r(l + (x/5) 2 ], where X = 5tan(0) is a stochastic 
variable representing the point where a foton released from S at the (stochastic) angle 8 crosses L. Now inserting n — 1 
lines, Li, parallel to L, between the light source, S, and L, we can apply Huygens principle, saying that Li will act as 
a line of point sources, where the light follows the path S — > X\ — > X\ + X2 — > ... — > X, where X = X\ + ... + X n . The 
variables Xi are independent and Cauchy distributed with scale parameter S/n. Thus Huygens principle is consistent 
with the fact that X = Q^iLi Xi)/n has the same distribution as Xi. 
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FIG. 7. Huygens principle applied to the light from a point source on a line. This illustrates the behavior of the averaging 
of Cauchy distributed stochastic variables. 



APPENDIX D: SIMULATIONS 



The simulations performed in this work only involves Cauchy noise, which is easily obtained from a random variable, 
X, uniformly distributed in the interval [— 7r/2, 7r/2], as Y — tan(X). More generally a random variable with an a- 
stable distribution |l3| is obtained from, 

Y = [sin(oA:)/cos(A:) (1/a) ] x [- cos([l - a]X)/ \og(W)} (1 ~ a)/a (Dl) 

where X is defi ned as above and W is another random variable uniformly distributed on the interval [0,1]. When 
simulating ( |2.lD by a discrete numerical time stepping the (fixed size) time steps usually needs to be much smaller 
than would be expected from numerical integration of the drift term alone. This is due to the large excursions from 
the tails of the distribution of the noise. It is thus important to use a stable integration routine for the drift term. A 
simple durable routine, which is the one used in these simulations, is Heun's integration scheme. The simula tion is 



performed as x(t + At) = x{t) + (f[x(t)\ + f[x(t) + f[x(t)At])At/2 + aAtl/ar)(t), where r)(t) is generated by (Dl) 
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